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Abstract 

We consider a general formulation of gradient flow evolution for problems whose 
natural framework is the one of metric spaces. The applications we deal with are con- 
cerned with the evolution of capacitary measures with respect to the 7-convergence 
dissipation distance and with the evolution of domains in spectral optimization prob- 
lems. 
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1 Introduction 

Shape optimization problems received a particular attention from the mathematical com- 
munity in the last years, both for the several applications that require the design of efficient 
shapes (for instance in Structural Mechanics and Aerospace Engineering) and for the difficult 
mathematical problems that have to be solved in order to obtain the existence of optimal 
solutions. In a very general form, shape optimization problems can be written as minimum 
problems like 

mm{F(tt) : tteA} (1.1) 

where A is a suitable family of admissible domains and F is a suitable cost function denned 
on A. Problems of this kind arise in many fields, and we quote the recent books [H |3l IH [T71 
[T8| |2TI [22] . where the reader can find all the necessary details and references. 

It is well known that the existence of optimal shapes only occurs in very particular situ- 
ations, where either some quite severe geometrical constraints are imposed to the admissible 
domains of the class A (like for instance convexity), or where the cost functional satisfies 
suitable monotonicity conditions (as it happens in several spectral optimization problems). 
When the existence of optimal shape fails, one has to deal with relaxed solutions, that be- 
long to a space much larger than the one of classical domains, and describe efficiently the 
behaviour of minimizing sequences for problem (11. ip . 

In this paper we are interested in problems of the form ( 11.11) arising in spectral optimiza- 
tion: the admissible class A is made of domains of M. d and the cost functional F is of one of 
the following types. 

Integral junctionals. Given a right-hand side / we consider the PDE 

-Au = /infi, u G Hq{VL) 
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which provides, for every admissible domain Q C M. d , a unique solution uq that we assume 
extended by zero outside of Q. The cost F(Q) = J(uq) is obtained by taking 



for a suitable integrand j. 

Spectral junctionals. For every admissible domain Q we consider the Dirichlet Lapla- 
cian —A which, under mild conditions on Q, admits a compact resolvent and so a discrete 
spectrum X(Q). The cost is in this case of the form 



for a suitable function $. For instance, by taking $(A) = we may consider the optimiza- 
tion problem for the k-th eigenvalue of —A: 



We will summarize some known facts about the minimization problems above and we 
deal with the problem of studying the shape evolution Q(t), starting from a given domain 
Qq according to a suitable definition of gradient flow. The theory of gradient flows in metric 
spaces has been recently developed in a great generality (see [2]) and in some situations it 
can be easily adapted to our purposes, in particular when we deal with relaxed problems. 
The extra compactness of the latter problems is of great help for proving the existence of a 
relaxed flow. On the counterpart, the flow is made of relaxed domains (capacitary measures 
in our case) and not of classical domains; we will show some examples in which, even starting 
from a very smooth initial domain Qq, the gradient flow quits the original admissible class 
A to evolve in the class of relaxed shapes (see [61 HE]). 

However, for some particular cases of cost functionals F, a different gradient flow can be 
considered; this will be made in Section H] where we show that a careful use of monotonicity 
properties of F allows to obtain an evolution path fi(t) made of classical domains. Some 
properties of the path Q(t) are studied, and some interesting open problems are pointed out. 

Let us close this introduction by mentioning that the idea of considering shape flows 
is somewhat reminiscent of many classical numerical treatments of ( II. ip where an initial 
(tentative) shape Qq is iteratively improved towards minimization. A numerical gradient 
flow perspective has in particular already been considered in [151 I2H] in connection with 
some applications to image segmentation, optimal shape design, and surface diffusion. 

2 Preliminary tools 

2.1 Capacity and quasi-open sets 

In the following we use the well-known notion of capacity for a subset E of M d : 



where Ue is the set of all functions u of the Sobolev space H 1 (M. d ) such that u > 1 almost 
everywhere in a neighborhood of E. If a property P(x) holds for all x G E except for the 




F(fi) = $(A(fi)) 



mm 



{A fc (fi) : tteA}. 
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elements of a set Z C E with cap(Z) = 0, we say that P(x) holds quasi- everywhere (shortly 
q.e.) on E, whereas the expression almost everywhere (shortly a.e.) refers, as usual, to the 
Lebesgue measure. 

A subset Q of M. d is said to be quasi-open if for every e > there exists an open subset Q £ 
of M. d , such that cap(f2 e Af2) < e, where a denotes the symmetric difference of sets. Actually, 
in the definition above we can additionally require that Q C Q £ . Similarly, we define quasi- 
closed sets. The class of all quasi-open subsets of a given set D will be denoted by A(D). In 
the following we always consider subsets Q of a bounded open set D <zM. d . 

A function u : R d -> M is said to be quasi- continuous (resp. quasi-lower semicontinuous) 
if for every e > there exists a continuous (resp. lower semicontinuous) function 
such that cap({w ^ u e }) < s. It is well known (see for instance [23]) that every function 
u G i/ 1 (lR rf ) has a quasi-continuous representative u, which is uniquely defined up to a set 
of capacity zero, and given by 

u(x) = lim / u(y)dy. 

e^o \B{x,e)\ J B{X;£) 

In the following we always identify, by an abuse of notation, a Sobolev function u with its 
quasi-continuous representative u, so that a pointwise condition can be imposed on u(x) for 
quasi-every x. In this way, we have for every subset E of M. d 



cap(E) = min 




By the identification above, a set f2 C M is quasi-open if and only if there exists a function 
u G H\R d ) such that SI = {u > 0}. 

The definition of the Sobolev space Hq(Q) can be extended for a quasi-open set fl; it is 
the space of all functions u G H l (M. d ) such that u = q.e. on M d \ fi, with norm 

IMIi^(n) = IfIIh 1 ^)- 

Most of the well-known properties of Sobolev functions on open sets extend to quasi-open 
sets. In particular, for every / G L 2 (D) there exists a unique solution of the PDE formally 
written as 

- Aw = / in fi, u G H£(Q) (2.1) 

that we consider extended by zero on D \ Q. The precise meaning of the equation above has 
to be given in the weak form 

ueH%(Q), [ VuVvdx= [ fvdx Vi> G Hq(Q) 

J D JD 

which turns out to be equivalent to the minimization problem 

min | J (^|Vtf -fvjdx : v G ^o(^)}- 

We denote the unique solution u of the problem above by Hn(f), which defines in this way 
the resolvent operator TZq. 
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2.2 7 _conver g ence an d ^-convergence 

The class A(D) of all quasi-open subsets of D can be endowed with a convergence structure, 
called 7- convergence. 

Definition 2.1. We say that a sequence of quasi-open sets (f2 n ) in A(D) --/-converges to a 
quasi-open set Q £ A(D) if for every f £ L 2 (D) we have that TZn n (f) converge to TZn{f) 
weakly in Hq(D). 

The following facts for the 7-convergence can be shown (see for instance [1]). 

1. In Definition 12. II it is equivalent to require the weak Hq(D) convergence only for / = 1. 
In addition, the 7-convergence of Q n to Q is equivalent to the T-convergence (see [ID] ) 
of the functionals 

Vm| 2 dx if u £ H (Q n ), +00 otherwise 

to the functional 

Vu| 2 dx if u £ Hq(Q,), +00 otherwise 

with respect to the L 2 (D) topology. 

2. It can be proven (see [I]) that if Q n — > Q in the 7-convergence, the convergence of 
the resolvent operators Hn n to IZn is in fact in the C(L 2 (D)) operator norm. In 
particular, the spectrum of 7Zn„ converges (componentwise) to the spectrum of TZn, 
hence the spectrum of —A on i/g(f2„) converges (componentwise) to the spectrum of 
-A on H^(n). 

3. The 7-convergence is metrizable on A(D); an equivalent distance to the 7-convergence 
is given by 

d 7 (fi!,fi 2 ) = \\n ni (i)-n n2 (i)\\ L 2 {D) . 

The 7-convergence is not compact; indeed it is possible to construct a sequence (fl n ) of 
domains such that the corresponding solutions TZn n (l) do not converge to a function of the 
form TZq(1). The first example of such a sequence was provided by Cioranescu and Murat in 
[9J by removing from the set D a periodic array of balls of equal radius r n — > 0. If the radius 
is suitably chosen they proved that the weak Hq(D) limit of TZn n (l) satisfies the PDE 

-An + cu = 1 in D, u £ H](D) 

where c > is a constant, and thus the sequence of domains (fi n ) cannot 7-converge to 
any domain Q. This is why, in order to study the behaviour of of minimizing sequences of 
domains, a relaxation procedure is needed. 

The relaxed form of a Dirichlet problem like (12. ip . has been obtained by Dal Maso and 
Mosco in [IT] where it is proven that the compactification of the metric space (A(D), d 7 ) is 
the set Mq(D) of all nonnegative regular Borel measures /1 on D, possibly +00 valued, such 
that 

fi(B) = for every Borel set B C D with cap(-B) = 0. 



/ 



/ 
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Note that the measures fi G M.q{D) are not finite, and may take the value +00 on large 
parts of D. For instance the measure 



oo D \ n (E) 



if cap(£ \ Q) = 0, 
-00 otherwise 



belongs to the class M.q(D). 

Given fi G M.q(D) we consider the space X^(D) of all functions u G Hq(D) such that 
J D u 2 dji < 00, endowed with the Hilbert norm 



\nx.{D) 



^ J |Vw| 2 dx + J u 2 dfi 



1/2 



This allows us to consider the relaxed form of a Dirichlet problem, formally written as 

-Au + jj,u = f in D, u G X^(D) 
and whose precise meaning is given in the weak form 

ueX^D), / VuVvdx + / uvd/2= / f(x)vdx Vi> G X^D). 

J D J D JD 

By the usual Lax-Milgram method we obtain that, for every fi G Aio(D) and every / G 
L 2 (D), there exists a unique solution u = TZ^f) of the equation above, which defines the 
resolvent operator 1Z^. 

If Q G ^4(-D) and fi = ooo\n then the space X^D) coincides with the Sobolev space 
Hq(Q) and TZn(f) = ^(/)- If / > 0, then by maximum principle the solution TZ^if) is 
nonnegative too, and then also / + Au = fiu is nonnegative. On the other hand, if / > 
we can write fi = (f + Au)/u which gives fi once u is known; of course we have /1 = +00 
whenever u = 0. Therefore, working with the class M.q(D) is in this case equivalent to work 
with the class of functions {u G Hq(D), u > 0, Au + f > 0}, which is a closed convex subset 
of the Sobolev space H^(D). 

The 7-convergence can be extended to the relaxed space A4 (D): we have fi n — > fi in 
the 7-convergence if for every / G L 2 (D) (it is equivalent to require it only for / = 1) 
the solutions TZ^Xf) converge to 7?. M (/) weakly in Hq(D). The main properties of the 
7-convergence on the space Ai (D) are listed below. 

1. The space Mq(D) endowed with the 7-convergence is a compact metric space; an 
equivalent distance to the 7-convergence is 

d 7 Oi,/i 2 ) = \\TZ^(1) -TZ^{1)\\ L 2 {D) . 



2. The class A(D) is included in Aio(D) via the identification Q 1— > ooo\n and A{D) is 
dense in Aio(D) for the 7-convergence. Actually also the class of all smooth domains 
Q is dense in A^o(-D)- 

3. The measures of the form a(x) dx with a G L l (D) belong to Aio(D) and are dense 
in A4q(D) for the 7-convergence. Actually also the class of measures a(x) dx with a 
smooth is dense in M.q{D). 
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4. If n n — y n for the 7- convergence, then the spectrum of the compact resolvent op- 
erator lZ^ n converges to the spectrum of 1Z^; in other words, the eigenvalues of the 
Schrodinger-like operators — A + \x n defined on X^ n (D) converge to the corresponding 
eigenvalues of the operator — A + \i. 

Since the 7-convergence is not compact, in order to treat shape optimization problems 
it is useful to introduce (see jl]) a convergence much weaker than 7, that makes the class 
A(D) compact. We call weak 7 this new convergence and we denote it by wry. 

Definition 2.2. We say that a sequence (f2 n ) of domains in A{D) wry- converges to a domain 
Vl G A{D) if the solutions w^ n = 1Zn n (l) converge weakly in Hq(D) to a function w G H^(D) 
(that we may take quasi- continuous) such that Q = {w > 0}. 

We list below the main properties of the w;7-convergence on the space A(D); for all the 
related details we refer the reader to [4]. 

1. We stress the fact that, in general, the function w in Definition 12.21 does not coincide 
with the solution wq = 7^(1); this happens only if Q n 7-converges to Q, which in 
general does not occur, because 7-convergence is not compact on A(D). 

2. The wy- convergence is weaker than the 7-convergence and, by its definition, it is 
compact, since the sequence wn„ = TZn„(^) is bounded in Hq(D) so it always has a 
subsequence {Q nk ) weakly converging to some function w G H^(D). 

3. If / G L^i^D) is a nonnegative function, then the mapping Q y j n fdx is wry-lower 
semicontinuous on A(D). 

4. If F : A(D) — > [—00, +00] is a 7-lower semicontinuous shape functional which is 
monotone decreasing with respect to the set inclusion, then F is w^-lower semicon- 
tinuous. For instance, integral functionals like f D j(x,un)dx with j(x, •) decreasing, 
where uq = 7Zn(f) and / > 0, and spectral functionals like $(A(f2)) with $ increasing 
in each variable, are w^-lower semicontinuous. 

2.3 Minimizing movements 

We recall here some notions and results in the direction of variationally-driven evolutions in 
metric spaces. In particular, we shall mention generalized minimizing movements and curves 
of maximal slope and their relation with gradient flows whenever a Hilbertian structure is 
available. In particular, we summarize concepts and results of interest for our purposes, 
referring to [2] for further details. 

In all of the following, (X, d) is a complete metric space, Wo G X is an initial condition, 
and F : X — >] — 00, +00] is a proper functional defined on X with effective domain D(F) = 
{x G X : F(x) < +00}. We let r(d) denote the topology in X induced by the metric d. 

At first, we shall mention the so-called minimizing movements theory which was intro- 
duced by De Giorgi in [12] in order to study evolution problems with an underlying varia- 
tional structure. The framework of the theory is very general and applies both to quasistatic 
evolutions as well as to gradient flows, under rather mild assumptions. 
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For every fixed e > the implicit Euler scheme of time step e and initial condition uq 
consists in constructing a function u £ (t) = w([t/e}), where [•] stands for the integer part 
function, in the following way 

{cP(v win)} 
F(v) H 

Definition 2.3 (Minimizing movements). We say that u : [0, T] — > X is a minimizing 

movement associated to the functional F and the topology r, with initial condition Uq, and 
we write u G MM(F, r, u ), if 

u e (t)4u(f) Vte[0,T]. (2.2) 

If the latter convergence holds for a subsequence e n — > 0, we say that u : [0, T] — > X is a 
generalized minimizing movement and we write u G GMM(F, r, Uq) . 

We say that a trajectory u : [0, T] — > X belongs to AC p (0, T; U), p G [l,oo], if there 
exists m G L p (0,T) such that 

d(u(s),u(t)) < I m(r) dr for all < s < t < T. (2.3) 

J s 

One can prove that, for all u G AC p (0, T; X), the limit 

\u'\(t) = \un dHs)Mt)) 
1 lw \t-s\ 

exists for a. a. t G (0,T). This limit is usually referred to as the metric derivative of u at t. 
In particular, the map t h- >■ |it'|(t) turns out to belong to L p (0, T) and is minimal within the 
class of functions m G L p (0,T) fulfilling Q, see [21 Sec. 1.1]. Let us recall [21 Prop. 2.2.3, 
p. 45] the following. 

Theorem 2.4 (Existence of generalized minimizing movements). Let the sublevels of F be 
r-compact in X. Then, for every initial condition uq G D(F) the set GMM(F,t,uq) is 
non-empty. Moreover, we have that GMM(F,t,u ) C AC 2 (0, T; X). 

2.4 Curves of maximal slope 

We say that a function g : X — > [0, +oo] is a strong upper gradient for the functional F if, 
for every curve u G AC(0, T; X), the function g o u is Borel and [21 Def. 1.2.1, p. 27] 

|F(n(t)) - F(u(s))| < / g{u{r))\u'\{r) dr for all < s < t < T. (2.4) 

J s 

In particular, if g is a strong upper gradient for the functional F and g o u G L 1 (0,T) we 
have that F o u is absolutely continuous and 

\{Fou)'\ < (gou)\u'\ a.e. in(0,T). 



7 



Definition 2.5 (Curve of maximal slope). Let g : X — > [0, +00] be a strong upper gradient 
for F. A trajectory u G AC(0,T; X) is said to be a curve of maximal slope for F with 
respect to its strong upper gradient g if 

-(Fou)'{t) = \u'\ 2 {t)=g 2 {u{t)) for a.e. t G (0,T) (2.5) 

In particular, F ou is locally absolutely continuous in (0,T), gou G L 2 (0,T), and the energy 
identity 



\J \u'\\r)dr + \f s 



t 

2 



2 

holds in each interval [s,t] C (0, T). 



g\u(r)) dr + F(u(t)) = F(u(s)) (2.6) 



The notion of curve of maximal slope is the natural extension to metric spaces of gradient 
flows in the Hilbertian setting. In particular, in case X is a Hilbert space with scalar product 
(■, ■) and norm || • || and F is, say, Frechet different iable one can readily check that 

v! + DF(u) = <=^> - \\u' + DF(u) || 2 = 

\\\u'\\ 2 + l\\DF(u)\\ 2 + (DF(u) 1 u'} = 
^ _(F o u)' = \\u'f = \\DF{u)f. (2.7) 

Hence, in the case of a smooth functional F the two notions of gradient flow and curve of max- 
imal slope coincide. More generally, curves of maximal slope in a Hilbert space correspond 
to gradient flows whenever some mild assumption is made on the Frechet subdifferential OF 
of F. The latter subdifferential is defined at points u G D(F) as 

r,„, s , F(w)-(F(u) + (v,w-u)) 
v g dF(u) =^ limsup ^ > 

with D(dF) = {u G D(F) : dF(u) ^ 0}. In particular, we have the following [2J Prop. 
1.4.1, p. 34]. 

Proposition 2.6 (Curves of maximal slope = gradient flows). Let (X,d) be a Hilbert space 
endowed with its strong topology. Moreover, assume that dF(u) is weakly closed for every 
u G D(dF). Then, u is a curve of maximal slope for F with respect to v 1— > \\d°F(v)\\ if and 
only if 

1 1 — ^ F(u(i)) is a.e. equal to a non-decreasing function, 
u'(t) + d°F(u(t)) 3 for a.e. t G (0,T), 

where d°F(u) is the subset of elements of minimal norm in dF(u). 

In the general metric setting the natural candidate for serving as a strong upper gradient 
for F is the local slope (see [21 [HJ [13]) of F defined at u G D(F) as 

|9F|(„)= limsup «")-^)) + 



d(u, v) 
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Note that indeed the local slope plays the role of the norm of dF (see ( 12. 7p ). In particular, 
in case X is a Hilbert space and F is Frechet differentiable at u G D(F), then \dF\(u) = 
\\dF(u)\\. 

In general, the function u \- > \dF\(u) cannot be expected to be lower semicontinuous. On 
the other hand, semicontinuity is crucial in order to possibly pass to the limit in (12. 6p (or, 
rather, in its time-discrete analogue). A way out from this obstruction consists in restricting 
the analysis to X-geodesically convex functionals. In particular, we say that a trajectory 
7 : [0, 1] — > X is a constant- speed geodesic if 

d( 7 (s), j(t)) = (t- s)d(j(0), 7 (1)) V0 < s < t < T 

and that a functional F is X-geodesically convex for A 6 K if, for all Uq,Ui G D(F), there 
exists a constant-speed geodesic 7 with 7(0) = u and 7(1) = U\ such that 

F(j(t)) < (1 - t)F(u ) + tF( Ul ) - ^(1 - t)d 2 (u , Ul ) Vt G [0, 1]. 

In case X has a linear structure, we shall simply (and classically) refer to the latter convexity 
condition as X-convexity. 

If F is A-geodesically convex and r(c?)-lower semicontinuous then [21 Cor. 2.4.10, p. 54] 
the local slope \dF\ is a strong upper gradient for F and it is r(d)-lower semicontinuous as 
well. The same holds if we relax the geodesic convexity assumption above by asking for the 
weaker property 

for all t>o, v i ^ D(F) there exists a curve 7 with 7(0) = Vq and 7(1) = V\ such that 

v H- $(£, v ,v) := —d(v,v ) + F(v) is (e" 1 + A)-convex on 7 for all < e < 1/A~ (2.8) 

along with the convention 1/A~ = +00 for A > 0. 

In particular, this entails the following [2j Thm. 2.3.3, p. 46]. 

Theorem 2.7 (GMM are curves of maximal slope). Let F fulfill the convexity assumption 
(I2.8P being r(d) -lower semicontinuous, and coercive, namely 

3e* > 0, u* G X : inf $(e*, «*, •) > -00. 

Then, given u G D(F), every u G GMM(F,T(d),u ) is a curve of maximal slope for the 
functional F with respect to its strong upper gradient \dF\. 

A suitably strengthened version of the convexity assumption (12.81) provides the possibility 
of proving a generation result. In particular, we shall be dealing with the following 

for all I**, vo, v\ G D(F) there exists a curve 7 with 7(0) = v and 7(1) = v\ such that 
v 1 — y &(e,v*,v) is (e:" 1 + A)-convex on 7 for all < e < 1/A~. (2.9) 

Note that property ( 12.91) is stronger than the former ( 12. 8ft for the latter follows as a particular 
case with v* = vq. On the other hand, (12. 9p combines in a crucial way the geodesic convexity 
properties of the functional and of the curvature properties of the underlying metric space 
[2]. One can in particular check that (12. 9p ensues if, in addition to the A-geodesic convexity 
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of F one requires v — > d?(v,v*) to be convex for all v* £ D(F). This is indeed the case 
of geodesically convex functionals on non-positively curved metric spaces such as Hilbert 
spaces [IH]. Note that, in the setting of assumption (I2.9p . no compactness of the sublevels 
of F is needed in order to prove the existence of curves of maximal slope and we have the 
following [3 Thm. 4.0.4, p. 77]. 

Theorem 2.8 (Generation of the evolution semigroup). Let F fulfill the convexity assump- 
tion (12.91) being r(d) -lower semicontinuous and coercive. Then, for any given u £ D(F) 
there exists a unique u = S(uq) £ MM(F, r(d) , uq) . Moreover, u is a locally Lipschitz curve 
of maximal slope for F with respect to its strong upper gradient \dF\, u{t) £ D(\dF\) for all 
t £ (0, T), and the map 1 1— > S(uo)(t) is a X-contraction semigroup, namely 

d(S(u )(t),S(v )(t)) < e- xt d(u ,v ) \fu ,v £ D{F). 

3 Curves of maximal slope of capacitary measures 

In this section we consider the compact metric space (M.o(D),d^) of capacitary measures 
endowed with the distance rf 7 introduced in Section 12.21 Let F : Ai (D) — >] — oo, +oo] 
be a 7-lower semicontinuous (relaxed shape) functional. By Theorem 12.41 for every initial 
condition fi £ M. Q (D) with /x £ D(F) there exists fi £ GMM(F,T(d 7 ), // ) and the discrete 
implicit Euler scheme reads 

£ argmin {f(jjl) + (3.1) 

The main purpose of this section is to study some properties of the generalized minimizing 
movement fi(t) and to see when it happens to be a curve of maximal slope. 
There is a natural one-to-one map between Aio(D) and the convex set 

X = {w £ Hl{D) : w>0, 1 + Aw > 0} C L 2 (D), (3.2) 

given by 

^ /,\ • i • 1 + Aw 
fx i — y i = iZ^l), with inverse w i->- fj, w •= . 

Moreover, the metric structure on Ai (D) and X is the same, since 

dybli, fa) = -W^\\ L 2 {D y 

Therefore, every functional F : A4 (D) — >] — oo, +oo] can be identified with a functional 
J : L 2 (D) ->] - oo, +oo] with D(J) C X by 

F{jj) = J{w^) or, equivalently, J(w) = F(fi w ). 

The variational flow for F in Ai (D) can be then obtained through the gradient flow of J 
in L 2 (D), generated by the implicit Euler scheme 

< +1 £ argmin | JH + ^- J \w - w n £ \ 2 dx}. (3.3) 
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Theorem 3.1 (Monotonicity). Assume that J : X — >] — oo, +00] is decreasing, in the sense 

Wi,w 2 G D(J), wi < w 2 a.e. =>- J(wi) > J{w 2 ). (3.4) 

Then, every w G GMM{J,T{dt,2^) ,w Mo ) zs increasing, in the sense that 

ti < t 2 ==>■ w(ti) < w{t 2 ) a.e. 

Proof. Let J be monotone in the sense of (13.411 . Given u>™, in the incremental step f ]3. 31) we 
have to solve the minimum problem 



min I J(w) + — I 

^ ^ E J D 



\w — w™\\ 2 dx !>. 



For every w G X, the function ma,x{w,w™} still belongs to X, since the maximum of two 
subharmonic functions is also subharmonic. Relying on the monotonicity of J we then have 
that 

J(max{w, W™}) + ^- J || max{u?,u?™} - w™\\ 2 dx < J(w) + ^- ^ ||w - w"|| 2 c?x, 

the inequality being strict as soon as ma,x{w,w™} 7^ to. Consequently, any minimizer w : = 
tt)" +1 should satisfy w > a.e., and thus any discrete flow is increasing. Passing to the 
limit as the step size goes to zero, we obtain that any generalized minimizing movement is 
increasing. □ 

Remark 3.2. We underline that the monotonicity assumption ( 13. 4 p is not equivalent to the 
monotonicity of measures. If \i\ < fi 2 in the classical sense of measures or, weaker, in the 
sense 

/ <f 2 (x)dfa < [ ^ 2 (x)dfi 2 VpeH^D), (3.5) 

J D J D 

then > w M2 q.e. The converse is in general false; here is an example: 

[l\ = ooLb(o,i) c 5 ^2 = ^B(p,i)dx + OO [b{0,R) c , 

where R is large enough, such that w^ 2 > q.e. Clearly, fii ^ fa- 
Remark 3.3. A typical functional satisfying the monotonicity assumption is a functional 
depending on w, of the form 

J(w) = / j(x, w(x)) dx, 
Jd 

where j : D x R — > R is continuous and decreasing in the second variable. In particular, we 
may take j(x, w) = —w which leads to the energy of the system for the constant force / = 1. 

Remark 3.4. If fx\ < fi 2 in the classical sense of measures, or in the weaker sense 03.5 j) . then 
it is easy to see that Xk(^i) < ^kifa)- We do not know if this is still true under the weaker 
(see Remark 13 . 2 j) condition that > q.e. 
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Example 3.5. An interesting question is the following: if we consider the generalized mini- 
mizing movement associated to the energy functional J(w) = — f D w(x) dx and start from 
a quasi-open set, will the flow remain in the family of quasi-open sets? 

As we show in the example below, by considering the topology of 7-convergence and 
allowing relaxation in general this does not happen, at least at the discrete level. This 
kind of phenomenon was numerically observed in the framework of quasi-static debonding 
membranes [6], where the evolution takes place in the family of relaxed domains. 

Let D = 5(0,2), tt = 5(0,2) \ 55(0,1) C R 2 and J{w) = -J D w(x)dx. Let us first 
notice that the mapping 

w J e (w) = — I w dx + — I \w — w \ 2 dx (3.6) 
Jd 2 £ Jd 

is strictly convex. We will prove that the solution w minimizing the first incremental step 
is of relaxed form, independently of the size of e > 0. For this purpose, we will first show 
that w is radially symmetric. It is not clear that the class X in (13. 2 p is stable by Schwartz 
rearrangement, in spite of the fact that we can use the convexity of the mapping above. 
Indeed, since Wq(x) = u (\x\) is radially symmetric, if w is a solution of the incremental 
step, then any rotation w o R of w, is also a solution. Using the strict convexity of (13. 6p . we 
conclude that for any rotation R the equality w = w o R holds, so w is radially symmetric. 

Assume now by contradiction that w corresponds to a non-relaxed domain, i.e. to a 
radially symmetric open set. This means that w = wq, where Q is a union of open annuli, 
centered at the origin. 

Denoting by A(s, t) the annulus B(0, t) \ B(0, s), with s < t, it is easy to see that optimal 
domains can only be of the form Q s = B(0,s) U A(s,2), which provide the corresponding 
radial solutions 

if < |x| < s 

w s (x) = u s (\x\) = { 1( A U|2 , , n2 log(|x|/2) \ m 

; log( S /2) J ~ 1 1 " 

In order to prove that relaxation occurs, it is enough to show that it is more effective to 
relax on dB(0, 1) the Dirichlet condition. In particular, given 

II = E% |_0B(O,l)+Oo|j9B(O,2), 

we have that F(/j,) is lower than any value F(Q S ). The corresponding solution w reads 

(s + u (\x\) if0<|x|<l 

and hence relaxation occurs at the first discrete step as soon as we prove that 

J £ (w) < J £ (w s ) VO < s < 2. (3.7) 
Indeed, by defining f(r) = min{l, log(r/2)/log(l/2)}, we have that 

J £ (w) = — 2ir / w (r)r dr — 2ixs i f(r)rdr + 7re / f 2 (r)rdr. 
Jo Jo Jo 
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Hence, relation (13 .7p is equivalent to 

/ (u s — u )r dr — — I \u s — u \ 2 r dr < e / frdr — ^-j f 2 rdr. (3.8) 
Jo 2e J Jq 2 J 

One can numerically check that the integral in the left-hand side above is non-positive (and 
indeed vanishes for s = 1 only). Hence, relation (13. 8 p follows, as we have that 

2 (f-\Ardr>0. 

Note that the latter argument is independent of e. As such, relaxation is expected to happen 
instantaneously for any generalized minimizing movement starting from Q . 

We shall assume that J : L 2 (D) — >] — oo, +oo] is A-convex, proper, and lower semicon- 
tinuous. For instance, J could be of the form 

J(w) = / j(x, w(x)) dx + Ix{w), (3.9) 
J D 

where j : D x R — > IR is a normal A-convex integrand and Ix is the indicator function of 
X, namely, Ix{w) = if w G X and Ix = +oo elsewhere. An example in this class is the 
torsional rigidity functional given by j(x,w) = w. Another example for J is 



J{w) = - I \Vw(x)\ 2 dx + / j(x,w(x))dx + Ix{w)- 
2 J d Jd 



Proposition 3.6. Let the Junctionals F : A4q(D) — >] — oo, +oo] and J : L 2 (D) — >] — oo, +oo] 
with D(J) C X be related by J(w) = F(fi w ) as above. We have the following 

a) F is X-geodesically convex if and only if J is X-convex, 

b) F fulfills (JZ1D if and only if J fulfills (12791) . 

Proof, a). Let F be A-geodesically convex. Given wo,w\ G D(J) there exists a constant- 
speed geodesic /i : [0, 1] — > Aio(D) such that = Wi, i = 0,1, and 

F(n(t)) < (1 - t)F(fi(0)) + tF(/i(l)) - ^(1 - t)rf 7 (/i(0), /i(l)). (3.10) 

By defining w(t) = w^), since 

\\w(s) - w(t)\\ L2(D) = <i 7 (/i(s),/i(t)) 

= (t-s)d 7 (/i(0),/i(l)) = {t-s)\\w -w 1 \\ LHD) V0 < s <t < 1, 

we readily have that w(t) = (1 — t)w + twi. By using (13. 10[) we conclude for the A-convexity 
of J as 

J(w(t)) = F(v(t)) < (1 - t)F(MO)) + tFQi(l)) ~\t{l- t)4(/i(0), /i(l)) 

A |2 



(1 - t)J( Wo ) + U( Wl ) - -t(l - t)\\ Wo - WiWbpy (3.11) 



13 



Assume now J to be A-convex, Fix /io, /ii £ D(F). By denning u>j = w ft for i = 0, 1 and 
iu(i) = (1 — t)u> + ^i we get that /i(t) = Hw(t) is a constant-speed geodesic joining /i and 
/ii as 

dy{n(s) t fi(t)) = \\w(s) - w{t)\\ L 2 {D) 

= (t - s)\\w - Wi||l2(23) = (t- s)d 7 (yU ,/ii) V0 < S < t < 1. 

Hence, by arguing exactly as in (13. lip the A-geodesic convexity of F follows. 

b). Let F fulfill (12. 9p and iw*, u>o, iui £ -D(J) be given. Define fi* = //„,,, /i = Hw , and 
A*i — A*toi; an d exploit (I2.9P in order to find the curve t /i(t) (possibly not a geodesic) 
joining /i and \i\ such that 

l^( M (t),^) < i^d^(0),^) + (1 -t)F(/i ) + ^d^^+tF^) 

-^^t(l-t)<(//o,/^i). (3.12) 

By letting iu(i) = u^t) we have that — u>*|| = <i 7 (/i(t), //*) and ||u>(£) — = 

d y (fj,(t), /i(s)) for all s,t £ [0,T]. Hence, relation (I3.12p entails that J fulfills (I2.9P as well. 

On the contrary assume that J fulfills ( 12.91) and that /i*, /io, /ii G D(F) are given. Define 
if* = wq = io Mo wi = and let t (-)■ w(t) be the curve whose existence is ensured by 
(12.91) . Then, by letting /i(i) = fi w (t) and arguing exactly as above we conclude that F fulfills 
(I2T9D as well. □ 

Proposition 13 .61 is based of the fact that the geometry of M.q and L 2 (D) coincide. Indeed, 
M o is a non-positively curved metric space. As such A-geodesic convexity in Ai implies the 
convexity property (I2.9P . We shall give a direct proof of this fact in the following. 

Proposition 3.7 (Geodesic convexity =>- (12.91) ). If F : Ai (D) — >]— oo, +oo] is X-geodesically 
convex then it fulfills the convexity property (12.91) 

Proof. Let /i*, /i , Hi G D(F) be given and define w* = w^, w = w^ , W\ = w(t) = 
(1 — t)w + twi, and /i(t) = Hw(t)- As /i is a constant speed geodesic we have that 

1 1 

//„//(*)) = — c/ 7 (/i(t),/iA + F(/i(t)) = — - w*||| 2(D) +F(fi(t)) 

+ (1 - t)F(fi ) + tF(/ix) - ^(1 - t)d 2 Jfi , //i) 
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1 + e\ 



[1 - *)$(e, /i*, /i ) + t$(e, /i*, /ii) - - ' ^ - t(l - t)c^(jio, Hi) 

whence the assertion follows. □ 

We shall now come to the existence results for evolution. Again, this can be formulated 
equivalently for trajectories of capacitary measures t i— > fi(t) £ Ai (D) or of their function 
representatives 1 1— > w(t) £ X. Let us start from measures. Theorem 12.81 yields the following. 
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Theorem 3.8 (Curves of maximal slope of capacitary measures). Let F : Ai (D) — >] — 

oo, +00] be proper, d 1 -lower semicontinuous and X-geodesically convex. Then, for any given 
fio G D(F) there exists a unique fi = S(fio) G MM(F, r(<i 7 ), /io). Moreover, fi is locally 
Lipschitz curve of maximal slope for F with respect to \dF\, u{t) G D(\dF\) for allt G (0,T) 7 
and 

d^S(/2 )(t), S(u )(t)) < dybto, v )e- xt V^o, ^0 e D[F) 

As for the function representatives w, the situation is that of classical gradient flows in 
Hilbert spaces. 

Theorem 3.9 (Gradient flow of subharmonic representatives). Let J : L 2 (D) — >] — 00, +00] 
with D(J) C X be proper, lower semicontinuous and X-convex. Then, for any given wq G 
-D(J) there exists a unique w = S(wo) G MM(F, r(<i i 2 (£>)), w ). Moreover, w is the unique 
solution of gradient flow of the functional J in L 2 (D). In particular, 

w' + dJ(w)30 a.e.in(0,T), w(0) = w . (3.13) 

Eventually, we have that 

\\S(w )(t)-S(v )(t)\\<e- xt \\w -v \\ \/w ,v eD{J). 

Note that, in case J is defined as in ( 13. 9 p via a smooth j, the inclusion in (I3.13P reads 

w' + d w j(x, w) + dl x {w) 3 

which is equivalent to w(t) G X and 

(w'(t) + d w j(-,w(t)))(w(t) -w)dx < a.e. in(0,T), Vu> G X. 

D 

More generally, in case J = M + Ix where M : L 2 (D) — y] — 00, +00] is a proper, convex, 
and lower semicontinuous functional with int D(M) n X ^ 0, the inclusion in (I3.13P means 
w(t) eXnD(M) and 



w'(t)(w(t) -w)dx < M(w) -M(w(t)) a.e. in(0,T), \/w E X f]D(M). 
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4 Variational flows of shapes: spectral optimization 
problems 

4.1 General shape evolution 

In this section, we deal with flows of shapes. There are no "standard" distances on the class 
A(D) of quasi-open subsets of D, and several choices will be studied in the sequel. A first 
natural distance is given by the Lebesgue measure of the symmetric difference set 
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Since two quasi-open sets may differ for a negligible set (think for instance in M 2 to a disk 
and a disk minus a segment), this is not a proper metric in A(D), so that one should consider 
equivalence classes in the family of shapes. 

For this purpose, for every measurable set M C D, we define the Sobolev space 

#o(M) := {u G H]{D), u = a.e. on D \ M}. 

We notice that, for a given open set Q, the space i^o(fi) may not coincide with the usual 
Sobolev space Hq(Q), the latter being possibly smaller, as soon as Q is non smooth. Nev- 
ertheless, for every measurable set M C D there exists a unique (up to a zero capacity set) 
quasi-open set Q(M) G .A(-D) such that 

^ 1 (M)=i7 1 (fi(M)). 

If Mi C M 2 then ft (Mi) C ft(M 2 ). In particular, ft(M) C M. Consequently, one can define 
the resolvent of the Laplace operator with Dirichlet boundary conditions associated to M 
by setting 

^ M : L 2 (D) -)• L 2 (D), ft A/ := ft n( M), 
and one can extend the 7 distance and the -^-convergences to measurable sets, by setting 

d 7 (Mi,M 2 ) = d 7 (ft(Mi),ft(M 2 )), M n ^ M if ft(M„) ^4 ft(M). 

Working with measurable sets instead of quasi-open sets in shape optimization problems 
associated to Sobolev spaces may, in general, severely change the result. Nevertheless, as 
soon as the functional F satisfies some monotonicity assumption, the problems become, in a 
certain sense, equivalent. We refer the reader to [5] for more details between this equivalence. 

So let us denote 

M(D) — {M C D : M measurable}. 

Let F : A(D) — >] — 00, +00] be a w^-lower semicontinuous functional, monotone decreasing 
for set inclusion and consider its extension to a functional on measurable sets given by 
F : M(D) ->] - 00, +00] 

F(M) = F(tt(M)). 

Functionals of the form 

F(fi) = $(Ai(n) J ...,A Jfe (fi)), 

where $ : M. k — >] — 00, +00] is increasing in each variable and lower semicontinuous are 
admissible, since for every k G N the fc-th eigenvalue of the Dirichlet Laplacian A&(f2) is 
decreasing with respect to the inclusion of measurable sets. As a consequence, F is monotone 
decreasing with respect to set inclusion as well. 

In M(D) the distance d c h ar is not compact. Nevertheless, relying on the monotonicity 
assumption on F, we have the following result. 

Theorem 4.1 (Generalized minimizing movements of shapes). Let F : A(D) — >] — 00, +00] 
be ivy-lower semicontinuous , monotone decreasing for set inclusion, and let M Q G M(D) 
such that Q(Mq) G D(F). Then, the set GMM(F,t w ^, Mo) is non-empty. Moreover, every 
M G GMM(F, t w1 , Mq) is increasing in the sense of set inclusion. 
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Proof. We prove first that every solution of the incremental Euler scheme 

M e " +1 G argmin A/eM(D) [f(M) + 1|M;aM| 2 }. (4.1) 



is increasing in the sense of inclusions. As a consequence of the monotonicity of F, we have 
that 

F{M U M £ ) + — |M"a(MU M £ n )| 2 < F(M) + — |M"aM| 2 , 

the inequality being strict as soon as \M £ \ M| is not zero. Consequently, every solution M 
of the incremental problem satisfies, if it exists, M £ C M. 

Let 0;^ = Q(Mk) Q where {M k )k is a minimizing sequence of measurable sets, each one 
containing M £ . By the compactness of the w~y- convergence, up to extracting a subsequence 
we have that uj k 10. As F(u) < liminf^oo F(w fc ), the measurable set wUM" is a solution 
to the incremental problem (I4.ip . since 

F(u U M k ) + U M fc )AM £ "| 2 < + U M,)aM;| 2 

= F(u) + i(wU M fc )AM £ "| 2 < liminf ( F(w fc ) + ^-|(o; fc U A4)aM;| 2 



< liminf f F(M fe ) + — \M k AM™\ 2 ) 
k-^+oo \ 2e J 



Finally, we set M £ n+1 :=uU M™. 

We rewrite the discrete flows in terms of quasi-opens piecewise constant sets t i— >■ u; e (i) 
and we pass to the limit as e — > 0. We reproduce in this setting the argument of [TJJ 
Thm. 3.2]. In particular, by using the monotonicity of the flows we have that the functions 
S s (t) = |o; e (t)Aa;o| are non-decreasing. Hence, by the classical Helly principle, at least for 
some not relabeled subsequence we have that S E (t) —> S(t) for all t G [0, +oo) and some 
non-decreasing function 5. The function 5 is continuous with the exception of at most a 
countable set of points N. We shall introduce the countable set M C [0, +oo) in such a way 
that 

G M, M is dense in [0, +oo), N C M. 

By a diagonal estraction argument and the compactness of the u>7-topology (still not re- 
labeling) one can find that oo £ (t) -? oj(t) for all t G M. Let us now fix t G [0, +oo) \ M, 
let t n G M such that t n /* t, and define uj(t) = U n u(t n ). We shall show that indeed u(t) 
coincides with the W7-Iimit of u) e (t). To this aim we exploit again the compactness of the 
u>7-topology, and extract a further (possibly t dependent) subsequence in such a way that 
oj e t(t) -4 u;*. By using the lower semicontinuity of the Lebesgue measure with respect to the 
u>7-topology [U Prop. 5.6.3, p. 125] We have that 



\u(t)Auj*\ = lim \u(t n )AUJ*\ < lim liminf \uj £ t (t n )Auj £ t (t)\ 

n— >+oo n— >+oo k— >+oo k k 



= lim liminf (S E t (t) - S e t (t n )) = lim (5(t) - S(t n )) 0. 

n— »+oo k— >+oo k k ' n _^_)_ 00 

Hence u(t) = a;*. In particular, the whole sequence u £ (t) u>7-converges to u(t) even for t ^ 
M. Finally, we have proved that t H- u(t) belongs to GMM(F, t w1 , Q(M )). Correspondingly, 
M £ has a pointwise limit M. In particular, M belongs to GMM(F,t w1 , Mq). □ 
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Example 4.2. Evolution of a ball. Let Qq = -8(0, Rq). For every e > 0, the discrete move- 
ment associated to Ai consists on balls. This is a consequence of the Schwartz rearrangement 
procedure. Consequently, the minimizing movement consists of an increasing evolution of 
concentric balls. 

For every given R > consider the ball -8(0, r), with r > R, which minimizes 

Al(B(0 , O) + 5«^ = r -A I (B(0, 1 )) + *^5!, 
where Ud denotes the Lebesgue measure of the unit ball in M. d . We obtain 

1^(5(0,1)) + ^ 
r 6 e 

which gives, for e small, 



-^Ai(B(0, 1)) + — {r d - R d )dr d - 1 = 



2A l( 5(0,l)) 
<Pu 2 d R 2d + l 

The radius R(t) during the evolution then satisfies the differential equation 

2A 1 (£(0,1)) 
{) d?u 2 d R? d + l 

which has the solution 

«(')-(c' t 4|j+ y"' 1,IM " 

Remark 4.3. F might be discontinuous. The functional F may be discontinuous on 
the curve t h- >■ Indeed, let us consider -F(fi) = Ai(O) and a generalized minimizing 

movement starting from Qq = -8(0, 1) U A^Ri^R?) where A(Ri,R<2) stands for the annulus 
centered at of radii 1 < R± < i?2- We choose Ri and i?2 in such a way that Xi(A(Ri, R2)) > 
B(0, 1). Hence, as the connected component A(R\, R2) of f2 does not contribute to the value 
F(Qq) of the functional, the intuition hint is that the generalized minimizing movement from 
flo will be Q(t) = A(Ri, Rq) U -6(0, fit)) with / increasing and discontinuous at f = R\. 

An interesting question is whether the evolution is stable in some particular classes of 
shapes. In particular, the interest in stability is related to compactness. In two dimensions 
of space, the class of simply connected open sets is compact with respect to 7-convergence. 
Moreover, in any dimension of the space, the class of convex sets is also compact with respect 
to 7-convergence. In the general case, we shall however remark that stability is not to be 
expected, as we argue below. 

Remark 4.4. Topological genus is not conserved. Let fi — ^(-^1,-^2) \C where C is a 
radial cut. As Ai(fio) > ^i(A(Ri, R2)), any generalized minimizing movement starting from 
Qo will immediately fill-in the cut so that the simply connected Qo gets to be non-simply 
connected. 

An example of an evolution from two simply-connected components to one non-simply- 
connected component is that starting from Q = B(0, 1)UU, where U is a suitable set, disjoint 
from 5(0,1). Let s = inf{r > 1 : 5(0, r) n U ^ 0} and assume that X^U) > Ai(5(0,s)). 
Since U does not contribute to Xi(B(0, r) U U) up to r = s, intuitively any generalized 
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minimizing movement starting from Q is of the form B(0, f(t)) U U with / increasing up to 
some intersection time. For suitably chosen sets U, after the intersection time the new set 
will not be simply connected. 

Note nonetheless that, due to the monotonicity of the flow, the number of connected 
components is non- increasing during the evolution. 

Remark 4.5. Convex shapes are unstable. Assume Qq to be the square [0,7r] 2 C M 2 . 
Consider the mapping T t : x G M 2 H- x + tv(x)n where v is suitably smooth and n is the 
outward unit normal to 8Qq and consider Q t := T t (Jlo). For the sake of definiteness, we shall 
normalize f dQo \v\ds = 1. We have that [TBI Thm. 5.7.1, p. 209] 







1 

fan 


dn 



v ds 



where U\ is the first eigenfunction (with unit L 2 norm) of the Dirichlet Laplacian on fl . 
Hence, by fixing a time step e > we can readily compute that the minimum of 



t^G v (t) -.= F(n t ) + hn t An \ 

le 





dui 
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fdflo 


dn 



v ds 



Ai(O ) + £ ( f \v\di 
Z£ \Jdn 



is attained at 



and corresponds to the value 
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v ds 



<dn 



dn 



v ds 



Ai(fi c 



In particular, as we readily compute that ui(x, y) = sin x sin y, the latter entails that in 
order to minimize v h- > G v (t v ) one would rather have v > and concentrate the mass of v in 
the middle of the sides of the square. Hence, we conclude that bumps are likely to develop 
from midpoints of the sides of the square so that convexity will be lost. 

We can formulate some open questions: 



Assume Qq is convex. Under which conditions on Qq an d m the minimizers of 

min{Ai(f2) : tto C tt, \Q\ — m} 



(4.2) 



are convex? According to the intuitive argument above, if Qq is a square and m is 
slightly larger than \£Iq\, then the optimal domains should not be convex. In [7] it is 
proved that for a "thin" rectangle Qq of sizes e and 1, and m < 7r/4, the solution of the 
shape optimization problem eqrefminxx cannot be convex, provided e is small enough. 

Let Qo be a convex set and assume that for every m > |Oo| every solution of the shape 
optimization problem ( 14.21) is convex. Is it true that then Qq is a ball? 



Is it true that the generalized minimizing movement associated to Ai in the framework 
of Theorem 14.11 will converge to a ball (rescaling if necessary)? 
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• Prove or disprove that the metric derivative of Ai computed at a bounded smooth set 
Q is given by 

oil 

Precisely, prove that 

Ai(fi)-Ai(n„) . dux 2 
hm sup < max — . 

Constraint on the measure. An alternative evolution, which does not require any rescal- 
ing, is to work in the class of sets with prescribed measure. Let c > 0. We consider only 
measurable sets M G M(D) such that \M\ = c. The incremental problem is given by: 

M £ " +1 G argmin MeM(D)i|A/|=c [f(M) + 1|M>M| 2 }, (4.3) 

No monotonicity can occur in this case, unless the flow is constant. The existence of a 
generalized minimizing movement associated to the incremental step ( 14. 31) is not clear. Nev- 
ertheless, one can construct discrete solutions of the incremental scheme. 

Indeed, let (Mk)k be a minimizing sequence in ( 14. 3p . We associate the quasi open sets 
Uk = O(Mfc) and, up to a subsequence, we can assume that uj n -^4 u;. If |u;| = c, then 1^ 
converges in L X {D) to l w , so that u is a minimizer. 

If \oj\ < c, we replace u with wU(7, where U is chosen such that \u U U\ = c and 

|M e n AM fc | -> |M £ n A(w UU)\. 

Consequently, u U C7 is a solution to the incremental step f 14 . 3 j) . 

An alternative way is to replace the measure constraint by adding a penalized term in 
the functional, i.e. to replace F(M) by F(M) + \M\. In this case, the existence of a solution 
to the incremental step relies on the lower semicontinuity of the Lebesgue measure for the 
^-convergence. 

Perimeter penalization. One can alternatively introduce a penalization on the perimeter. 
In this case, the incremental step reads 

M £ n+1 G argmin MeM(D) [f{M) + P D {M) + 1|M;aM| 2 }, (4.4) 

The topology given by d c h ar turns out to be compact on the sub levels of F + Pd- 

Hausdorff distance. There are several other geometric distances in the family of open 
sets, but they are hardly compatible with the 7-convergence. Nevertheless some partial 
observations can be done. 

Let due denote the Hausdorff complementary distance in the family of open subsets of 
D, given by 

djjc(rii, Q 2 ) — max \d(x, D \ Qi) — d(x, D \ f2 2 )|. 

Assume F is increasing with respect to the set inclusion. Then, there exists a solution of 
the iteration step 

mini' iZ H , 



dui 
dn 
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which is of the form 

n = D\(n c + B h ). 

Let fl C D be open. Then, f2flf2o is open and, as F is increasing with respect to set inclusion 
and d#c(f2 D VLq, Qq) < d^iP-, ^o), we have 

f(q n o ) + ^-4c(fi n o , ^o) < F(n) + — d^ c (n, ft )- 

Now, define h = d 2 Hc (Q D S7 ? ^o) an d observe that -D \ (fig + B^) C Q n f2 and dH^{D \ 
(n c + B h ),n ) = h. Then, 

F(D \ (fig + + \ («S + n ) 

< f(q n n ) + 7^d 2 HC (n n n , fio) < F(n) + — 4- c (fi, n ) vacd open. 

Finally, there exists a generalized minimizing movement associated to F and d^c of the 
form t ^ D\ (Qq + -B/(t)), where / is continuous and increasing. 

The same argument for decreasing functionals F and the Hausdorff distance 

d H (F u F 2 ) = max \d(x, F t ) - d(x, F 2 )\ 

can be repeated. The only point which is more delicate is concerned with the fact that the 
Hausdorff distance is not a "proper" metric in the family of open sets. Nevertheless, one 
can prove the existence of generalized minimizing movement associated to F and du of the 
form t i — y int(f2o + -B/(t))j where F is continuous and increasing. This solution relies on 
the equivalence relation in the family of the open sets: fii = Q 2 if &i = ^2 and on the 
redefinition of the Sobolev space 

H*(Q) := {u G Hq(D) : u = a.e. on D\U). 
4.2 Flows of convex shapes 

In this section we deal with the evolution of convex open sets. We introduce the family 

/C(-D) = {K C D : K open and convex}. 
There are different possible distances on K which have the same convergent sequences 

• the Hausdorff distance; 

• d 2 (Ki,K 2 ) = H&RT! — bK 2 \\L 2 (D)i where b K is the oriented distance function, defined by 
bx{x) = —d(x, dK) for x e K and b K (x) = d(x, dK) for x E D \ K; 

• the L 1 distance of the characteristic functions d c h ar (Ki, K 2 ) = f D \ l Kl — 1 K2 \ dx. 

A slightly different distance, defined on the equivalence classes of homotopic convex sets 
is the Fraenkel relative asymmetry, defined by 



A(K U K 2) := inf + **»)! I where A :- 

i GK n I \Ki\ J 





l/n 


\K 2 


1/n 
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Assume that F : K,(D) — > R is a 7-lower semicontinuous shape functional which satisfies 
F(K n ) — y +00 as soon as K n converges to a degenerate set. Since in the class of convex 
sets, the u^-convergence coincides with the 7-convergence, it is useless to require W7- lower 
semicontinuity. Notice that all previous topologies are compact on sublevels of F. By 
applying Theorem 12.41 we have the following. 

Theorem 4.6 (Generalized minimizing movements of convex shapes). For d = dH,d c h ar , 
or di, and for every initial convex set K 6 D(F), we have that GMM(F,r(d), K ) is non- 
empty. 
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